Tunneling into strongly biased Tomonaga-Luttinger liquid 
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We calculate the tunneling density of states for a Tomonaga-Luttinger liquid placed under a 
strong bias voltage. For the tunneling through a side-coupled point contact, one can observe the 
power law singularities in the tunneling density of states separately for the right- and left-movers 
despite the point-like tunnel contact. Deviations of the nonequilibrium tunneling exponents from 
the equilibrium case are discussed. 



Introduction. Power-law suppression of the tunneling 
density of states into the Tomonaga-Luttinger liquid is 
one of the most profound manifestations of interactions 
in one dimensional (ID) electron systems It paves 
the way to get the information about interactions in ID 
systems experimentally [l, Q. Exact calculation of the 
tunneling density of states is possible in the framework 
of the Tomonaga-Luttinger liquid (TLL) model using the 
bosonization method 0,13] ■ The basis of the TLL model 
is the linearization of one particle spectrum around the 
right and left Fermi points. At the same time, a proper 
account for the finite curvature of the one particle disper- 
sion is necessary for theoretical description of high energy 
excitations in ID electron systems and in particular for 
the description of strongly nonequilibrium ID systems 



In this letter we consider tunneling from the Fermi- 
liquid reservoir into the nonequilibrium TLL through a 
point tunnel contact (see Fig. [1]). The nonequilibrium 
conditions are created by a strong transport voltage V^d 
applied to a TLL channel. In the equilibrium the whole 
system is filled by electrons up to the Fermi energy Ep. 
Finite source-drain voltage results in the shift of the 
chemical potentials for the right- and left-movers to the 
quasi-Fermi energies Ep-\-eVsd/'2 and Ep — eVsci/'i respec- 
tively. At strong enough voltages, the nonlinearity of the 
electronic dispersion leads to different Fermi velocities of 
the right- and left-movers vu.l — \J (2i?F ± eV^d)/w*, 
as depicted in Fig. [T] (Here, m* is the effective electron 
mass.) In turn, the tunneling densities of states for the 
left- and right-moving spectral branches differ. Further- 
more, since the direction of partial tunneling currents 
into the left branch and out of the right branch are op- 
posite, these two tunnel currents do not compensate any 
more even at zero voltage V^c at the point contact (see 
Fig- mi- Therefore, a finite tunnel current flows between 
the nonequilibrium TLL and the reservoir. This current 
depends as a power law both on the source drain volt- 
age in TLL and on the voltage on the point contact 
Vpc, with the exponent reflecting the interaction strength 
in TLL. The exponent differs from the one describing 
the tunneling anomaly in the equilibrium TLL. We pro- 
pose a modification of the bosonization approach that 
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FIG. 1: (a) Occupation of the dispersion parabola in ID chan- 
nel in presence of the source-drain voltage T4d at negligible 
Vpc. The Fermi energy of the 2DEG lies nearly in the middle 
between quasi Fermi levels of the biased TLL. (b) Schematics 
of the device proposed. 



allows to calculate the nonequilibrium tunneling density 
of states in TLL analytically. The dependence of the cur- 
rent through the point contact on the source-drain and 
point contact voltages /pc(K5d,V^c) taking into account 
electron-electron interactions in the ID channel is the 
main result of this paper. 

To provide a theoretical description of TLL under a 
strong source-drain voltage, we have to deal with quite 
unusual ID electron liquid with different densities of 
states for opposite chiralities. To our knowledge, this 
modification of Tomonaga-Luttinger model has not been 
solved so far. The Hamiltonian for the TLL channel con- 
sists of a sum of kinetic and interaction energies 
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where the kinetic energy term takes into account different 
quasi-Fermi energies and different Fermi velocities for the 
left- and right-movers 
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dx 
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and Hmt is the standard charge density interaction 
Hamiltonian of TLL. In Eq. ^ the energy is counted 
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from the position of the Fermi energy in the unbiased 
TLL. For description of the tunnehng it proves conve- 
nient to perform a gauge transformation ijjR(L){x) ^ 
U {x)'^ ijj ji(^]^-^{x)U {x) . The matrix U can be represented 
as U{x) ~ {j ii{x)lJ l{x) with 



ixeVsd 



2hv 



Ur{l){x) = exp 

(3) 

This gauge transformation absorbs the quasi Fermi ener- 
gies into the spatial dependence of the transformed fields. 

Using the bosonization identity, the total Hamiltonian 
(fT|) can be written in a bosonized form often encountered 
in the literature, namely 
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where cj)^ = (0/j(a;), 0^ (a;)) is the bosonic field corre- 
sponding to the fermions in TLL. Here the terms with 
quadratic number operators N'^ ^ are omitted. The in- 
teraction constants can be found via Fourier transforms 
of a given interelectron potential For the spinless 
case considered here we have 174 = V(q — Q)/hvR and 
92 = [V{q = 0) - V{q = - k^)] /hvR, where V{q) is 
the Fourier transform of the screened interelectron inter- 
action potential, and kp'^ are the Fermi wave vectors. 
The constant g — vl/vr is assumed to be a measure of 
the dispersion law asymmetry. 

Diagonalization of the Hamiltonian. The conventional 
way to diagonalize the Hamiltonian of any interacting 
system in the framework of the Tomonaga-Luttinger 
model is to introduce so-called "dual fields" . The stan- 
dard couple of "dual fields" is defined as a difference 
and a sum between boson fields with opposite chiralities 
[^]. However, because of the different Fermi velocities for 
right- and left-moving electrons, the introduction of dual 
fields cannot be made in that direct manner for the prob- 
lem at hand. Rather, we introduce additional fictitious 
bosonic fields 4>'r{x) and in a way which does not 

change the dynamics of the system. Then, in terms of 
the four-component bosonic fields B = (/'l, f/)'^, 4''r)'^ i 
the Hamiltonian (|4]) can be rewritten as 
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(5) 

Note that there is no coupling between the original and 
fictitious fields, which guarantees that the dynamics of 
the original fields remains unchanged. Using the fact of 
the chiral symmetry between the original and fictitious 
branches we form the two couples of dual fields out of 
^_R,L and 



^'l r which read 



$1 
$2 
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Substituting ([6]), ([7]) into the Hamiltonian ^ we obtain 
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where 
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the upper and lower signs corresponding to Af$ and Mq 
respectively, and <i> = ($1, $2)^, B = (61, 62)'^. 

Hamiltonian ([5]) can be brought to the canonical form 
applying a composition of unitary rotations and rescal- 
ings that preserve conformal invariance. We first diag- 
onalize the matrix M$ = P$A$P^^ by a unitary rota- 
tion of the fields with a matrix P$. Then we rescale 
the fields while preserving duality relations (P$$)j — > 
(P**), /^Af, (P*e), ^ ^(P$e),, so that after the 
rescaling the fields $i are coupled by the unity matrix, 
and the Hamiltonian acquires the form 

H = ^J^ [(a.$^)i(5,a>) + (9,e^)Me(9.e)" , 

(10) 

where Me = %/A^P<j'^AfeP<E.%/iW- 

Further we diagonalize the quadratic form with the 
fields by a unitary rotation with a matrix Pq, i. e. 
Pg^AfePe = Ae. The unity matrix that couples the 
fields $ remains unaffected by this transformation. Fi- 
nally, we repeat the rescaling of the fields while preserv- 
ing duality and bring the Hamiltonian to the canonical 
form 
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where the new dual fields 77 — (771,772)"'", C ~ ("^1,^2)"^ are 
related to the original ones by transformations $ = S'$77, 
Q — SqS, with 



The new velocities Ui arc determined by Ui — vr^ 
and can be written explicitely in the form 
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Tunneling density of states and tunneling current. The 
tunneling current can be represented as a sum of par- 
tial currents between the 2DEG and the right- and left- 
moving branches of TLL 



^pcR 



^pcL 



(14) 
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with the partial currents given by [8[ 
47re , f 

Ipc, R(L] = — \ t\ de[nF{e- eVpc ± eVsd/2) 

-npie)] !^2DEG'^fl(L) (e) , (15) 

where t is the tunnehng matrix element describing the 
point contact, and rt^? is the Fermi distribution. The 
electron density of states in a 2DEG i^2DEG is just an 
energy- independent constant. The tunneling density of 
states i^fl(L)(e) is calculated from the one particle Green's 
function of an electron in the TLL. Using the bosoniza- 
tion identity Vacl) ^ e*'^«(^'), Eqs. ©, and the 

relations between the initial and canonical bosonic fields 
given by Eqs. (|12p . we obtain the imaginary time one 
particle Green's function in the form 

Gr{l){x,t;x,0) = ~{TiPr(l){x,t)tI}^^^^{x,0)) 
7-e75(*i(2)=Fei(2)),,,g-75(i>i(2)Tei(2))^„ 



= [(a/(ui|r|)]^-<-' [a/KIrl)]"-'-' , 



Performing analytical 



where = \ ^3%? + (5^)^^ 
continuation to real time and Fourier transformation [31 , 
we finally obtain the chiral tunneling density of states in 
the form 



^11(21) ^12(22) 
U-\ 



(16) 



sin [tT (6ii(21) + ^12(22))] T (l - &11(21) - &12(22)) , 

where a is a short-length cutoff, and lo ~ e/h. 

The nonequilibrium tunneling density of states can be 
best seen in the measurements of the differential conduc- 
tances (9/pc/9Vpc,sd at small voltages V^c on the point 
contact. The expression for the differential conductance 
can be obtained straightforwardly from Eqs. p4|) and 
(fTSll . The result at zero temperature can be written in 
the form 



i9Vpc,sd 

eV^d 



C Z '^2DEG X 



(17) 



where C = 1/2 for dlpc/dV^d and ^ = 1 for dIpc/dV-pc- 
In general, the dependence of the differential conductance 
on V"sd and on Fpc is smooth and is determined not only 
by the power-law singularity in l|16p but also by the de- 
pendence of g on Vsd- The latter makes the powers hij 
and the velocities wi,2 dependent on the voltage, in ac- 
cordance with Eq. (jlSp . The power-law singularities of 
nonequilibrium chiral densities of states can still be seen 
in the differential conductance at Vpc = ±T4d/2, as it 
follows from ([TT]). At these voltages, the Fermi level in 
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FIG. 2: (Color online) DiflFerential conductance dIpc/dVpc 
in units of Sttb^ !/2DEG'^iD/7r?i (i^id being the density of 
states at the Fermi level of the unbiased noninteracting ID 
channel) as a function of voltage at the point contact for dif- 
ferent values of 172 at a given 174. The bias voltage V^d is taken 
equal to O.lEp, and the other parameters are relevant for typ- 
ical GaAs-based electron gases: m * = 0.067me, a = 5.65A, 
Ef = 30meV. The conductance dependencies (?/pc/9Kid(Vpc) 
exhibit singularities at Vpc = ±14d/2, in accordance with 
Eq. ini). 



the Fermi liquid reservoir coincides with the quasi Fermi 
energy for the left- or the right-moving fermions in TLL, 
and the tunneling density of states in the corresponding 
channel is suppressed. These singularities are illustrated 
in Fig. [5] for dIpc/dVpc at different values of g2 which 
characterizes the screening of electron-electron interac- 
tions. At strong screening, when V{q = 0) ~ V{q = 
kp — kp) and hence 172 ^ 94, the differential conduc- 
tance exhibits a sharp dip close to Vpc = ±T4d/2. At 
weaker screening (172 closer to 54) the conductance is get- 
ting smaller, and the dip is essentially broadened. 

Furthermore, since 611 = 622 and 621 = 612 for any pa- 
rameters g and 172,4 the energy dependencies of un and 
vl are the same. There is however a difference in the 
prefactors that depend on the powers of plasmon veloci- 
ties ui and U2. This asymmetry is maximal at 612(21) = 
(i. e. when 172 = 0). In that case the only effect of in- 
teractions consists of the renormalization of the plasmon 
velocities, while the singularity in the tunneling density 
of states disappears. At larger g2 (i. e. 612(21) 7^ 0) 
the electrons with opposite chiralities interact with each 
other that leads to the alignment of the chiral tunnel- 
ing densities of states vr and i^l, as one can see from 
Eq. p6|) . The dependence of the tunneling exponent 
a — 611 -|- 612 — 1 = 621 + 622 — 1 on the asymmetry 
of the Fermi velocities is shown in Fig. [3l The power a 
increases with the grows of asymmetry (smaller g) from 
its equilibrium value at g = 1. Therefore, the singularity 
in the tunneling density of states becomes stronger with 
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FIG. 3: (Color online) The dependence of the tunneling ex- 
ponent a = foil + fei2 — 1 = &21 + 622 — 1 on the asymme- 
try of the Fermi velocities vr and vl expressed through the 
parameter g = vl/vr. The unbiased TLL corresponds to 
g = 1. The asymmetry increases [g diminishes) increasing 
the bias voltage Kd- The suppression of the tunneling density 
of states grows with bias voltage. Inset: Differential conduc- 
tance S/pc/c^Vpc as a function of voltage at the point contact 
for different values of Vsd at a given 32 = 54 = 0.5. The other 
parameters are taken the same as for Fig. [5] 



increase of the bias voltage Vsd ■ It is interesting that the 
finite size of the TLL also leads to increasing of the ex- 
ponent of a power in the tunneling density of states Q . 
Note, however, that it is much easier to change the bias 
voltage than the length of the wire. 

Another interesting feature of the system considered 
is the finite current through the point contact even at 
vanishing Vpc ■ This is again due to the chiral asymmetry 
of the plasmon spectrum subject to the bias voltage V^d- 
The existence of a finite current follows directly from 
Eq. (fT7)). but in oder to ease the understanding of that 
fact, we rewrite dl^d dVsd in the hmiting case of V^c = 
and strong screening when g2 = 0. The latter leads to 

bl2 = &21 = 0, bii = 622 = 1, Ml = Ufl(l -I- 34), U2 = 

Vl + viiQi, and the differential conductance assumes the 
form 



pc 



dVs 



sd 



\t\ v2n-E,G{vB.- Vl) 
h'^ {vL + VB.gi){vii + VRgi) ' 



(18) 



From Eq. (|T8|) it is clear that, on one hand, the tunnel- 
ing current into the TLL can be suppressed due to the 
electron-electron interactions. On the other hand, one 
can facilitate the tunneling applying the bias voltage V^d 
to the TLL. We emphasize that one does not need to 
change Vpc- Since high voltages at the point contact are 
not always possible in the linear response measurements, 
the biasing of the TLL might be a powerfuU tool to study 
its transport properties. 



Conclusions. In conclusion, we showed that in the ex- 
periment on the tunneling into a strongly biased TLL 
through a point contact, the power law singularities in 
the tunneling densities of states can be seen separately 
for the right- and left-movers. We obtained analytical ex- 
pressions for the chiral tunneling densities of states that 
turn out to be different from the equilibrium case. The 
predicted behavior can be observed in experiments with 
GaAs quantum wires of nominal width ^ 14 nm, where 
the application of a strong bias voltage does not cause the 
population of the next one-dimensional subband. Fab- 
rication of such wires lies within the range of current 
experiments [3, [13] ■ We also developed a method of di- 
agonalization of TLL Hamiltonian with different Fermi 
velocities for the right- and left-movers. The method can 
be usefuU for a number of problems which involve chiral 
asymmetry of the density of states such as a TLL wire 
in an external magnetic field. 
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